Normal scaling in globally conserved interface-controlled coarsening of fractal clusters 
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Globally conserved interface-controlled coarsening of fractal clusters exhibits dynamic scale invari- 
ance and normal scaling. This is demonstrated by a numerical solution of the Ginzburg-Landau 
equation with a global conservation law. The sharp-interface limit of this equation is volume pre- 
serving motion by mean curvature. The scaled form of the correlation function has a power-law tail 
accommodating the fractal initial condition. The coarsening length exhibits normal scaling with 
time. Finally, shrinking of the fractal clusters with time is observed. The difference between global 
and local conservation is discussed. 
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Dynamics of growth of order from disorder in coars- 
ening systems with long-range correlations is an intrigu- 
ing problem which appears in phase ordering and in 
many other apphcations. Systems with long-range corre- 
lations are often characterizable by fractal geometry [^,^ . 
Therefore, a lot of attention has been devoted recently to 
a variety of "fractal coarsening" problems [^|-^ . A typi- 
cal (though not the only) setting for fractal coarsening is 
the following. At an earlier stage of the dynamics a frac- 
tal cluster (FC) develops due to an instability of growth 
of the "minority phase" [§,0 . When the mass (or heat) 
source is depleted, coarsening by surface tension becomes 
dominant. The main question concerns possible scaling 
relations and universality classes in this type of coarsen- 
ing. 

A major simplifying assumption in the analysis of a 
coarsening process is dynamic scale invariance (DSI) (see 
review [|). The DSI hypothesis was first applied to frac- 
tal coarsening by Toyoki and Honda [Q who considered 
systems with non-conserved order parameter. Implica- 
tions of (mass) conservation in fractal coarsening were 
considered more recently |^,^ . Most remarkable of them 
is predicted shrinking of the FCs in the process of coars- 
ening. However, there has been no convincing evidence 
(neither in experiment, nor in simulations) in favor of 
DSI in conserved fractal coarsening. Moreover, anoma- 
lous scaling and breakdown of DSI were observed in re- 
cent simulations of locally conserved edge-diffusion- ||] 
and bulk-diffusion-controUed fractal coarsening. 

This communication reports our finding that DSI and 
normal scaling hold in the process of interface- controlled 
fractal coarsening with a globally conserved order param- 
eter. This system is apparently the first realistic con- 
served fractal coarsening system where this simplifying 
and beautiful concept is found to work. 

Globally conserved interface-controlled coarsening is 
accessible in experiment. Consider the sublima- 
tion/deposition dynamics of a solid and its vapor in a 
small closed vessel kept at a (constant) low temperature 
. As the acoustic time in the gas phase is short com- 



pared to the coarsening time, the gas pressure (and, con- 
sequently, density) remains uniform in space, changing 
only in time. This character of mass transport makes 
the coarsening dynamics conserved globally rather than 
locally, which leads to a different kinetics. Another ex- 
ample appears in the context of attachment / detachment- 
controlled nanoscale fluctuations at solid surfaces ||l8|,|l9| . 
There is also a strong recent evidence in favor of interface- 
controlled transport during the cluster coarsening in 
driven rapid granular flows [^o| . 

Here is an outline of the rest of the paper. We shall 
work with the Ginzburg-Landau equation with a global 
conservation law. The corresponding sharp interface the- 
ory is reducible, at large times, to volume preserving 
motion by mean curvature. Assuming DSI, one can 
then predict scaling behavior of the correlation function, 
shrinking of FCs with time and normal scaling of the 
coarsening length l{t). Our extensive numerical simula- 
tions of the coarsening of two-dimensional (2d) diffusion- 
limited aggregates (DLAs) support all these predictions. 
We shall conclude by pointing out the main difference 
between global and local conservation. 

We adopt a simple Landau free energy functional: 

F[u]^ J [{l/2)iVu)^ + V{u) + Hu] d'^r, (1) 

where V{u) = (1/4)(1 — u^)^ is a double-well poten- 
tial, u(r,t) is the order parameter and fluctuations are 
neglected. The effective "magnetic field" H = H{t) 
changes in time so as to impose a global conservation 
law: < u >= const, where < ... > denotes a spatial 
average: 

<. ..>=L-^ !(...) d'^r, (2) 

L is the system size and the integration is over the whole 
system. The dynamics is described by a simple gradient 
descent: 

^^J-£^^^u + u-u'-H{t). (3) 
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It follows from Eq. 



and the conservation law that 



Hit) 



--< u 



> (either no-flux, or periodic bound- 



ary conditions are assumed), so Eq. (|3|) is a nonlo- 
cal reaction-diffusion equation |21-24|. In the context 
of phase ordering it can be called the Ginzburg-Landau 
equation with a global conservation law. To make a the- 
oretical progress, one should work in the sharp-interface 
limit valid at late times, when the system already 
consists of large domains of "phase 1" and "phase 2" di- 
vided by a sharp interface. At this stage H{t) is both 
small, H{t) <g; 1, and slowly varying in time. The phase 
field in the phases 1 and 2 is uniform and rapidly adjusts 
to the current value of H{t), so u = — 1 — H{t)/2 and 
1 — H(t)/2, respectively. For brevity, we will consider 
2d-case. The normal velocity of the interface is 



(4) 



where s is the coordinate along the interface and k is 
the local curvature of the interface. The positive sign of 
Vn corresponds to the interface moving toward phase 2, 
while K is positive when the interface is convex towards 
phase 2. 

An equation for H{t) follows from the global conserva- 
tion law which takes the form 



4A(t) 
L2 



— H[t) — const , 



(5) 



where A{t) = ju{Tt)>o'^^^ cluster area. 

Eqs. (^) and (||) make a closed set and provide a gen- 
eral sharp-interface formulation to the Ginzburg-Landau 
equation with a global conservation law. Often they can 
be simplified further. Compute the area loss rate of the 
cluster: 



A{t) 



Vn{s, t)ds — K{t) 



—H{t)~K{s,t) 



, (6) 



where K{s,t) is the curvature averaged over the whole 
interface: 



K{s,t) 



1 



A(t) 



k(s, t)ds , 



(7) 



and A{t) is the cluster perimeter. If the cluster area is 
conserved, then H{t) — (\/2/3)k(s, i) which yields 



Vnis,t) = K(s,t) - K(s,t) . 



(8) 



This is area-preserving motion by curvature (in 2d), or 
volume-preserving motion by mean curvature (in 3d) 
p2| , p4|j25| . Dynamics (P) shortens the interface length 
(in 2d), or area (in 3d) p^ , p5[ |. This nonlocal coarsening 
model is simpler than the better known "Laplacian coars- 
ening model" (derivable from the Cahn-Hilliard equation 
[0,^) which describes the late-time asymptotics of the 
locally-conserved bulk-diffusion-controlled coarsening. 



Assuming DSI and using Eq. (^, one obtains normal 
scaling for the characteristic coarsening length: l{t) ~ 
t^/2 Therefore, global conservation does not change the 
scaling law. The same result (again, when assuming DSI) 
follows from dynamic renormalization group arguments 
applied to Eq. (^) (with a Gaussian white noise term) 
prf . For short-range correlations this result was sup- 
ported by particle simulations of critical and off- 
critical quench, and by a numerical solution of Eq. 
(^) for both critical, and off-critical quench |^ . 

Let us return to fractal coarsening. The initial con- 
ditions represent FCs that are characterizable by their 
fractal dimension D on an interval of scales between the 
lower cutoff Iq and the upper cutoff Lq. The DSI-based 
coarsening scenario assumes that the fractal di- 

mension of the cluster remains constant on a shrinking 
interval of distances between the lower cutoff l{t) (which 
has the same dynamic scaling as the coarsening length), 
and the upper cutoff L{t). Now, the perimeter A and 
area A of the EC can be estimated as H 



A-/(L/0^ and Ar^l^{L/l) 



D 



(9) 



respectively. Then area conservation yields L ^ 
l(D-2)/D _ ^-{2-D)/2D ^j^j^j^ imphes that the FC shrinks 
with time [|§. This follows K{t) ~ l-^{t) - i'^/^. 
One can also predict the asymptotic shape of the equal- 
time pair correlation function at large times: C{r,t) 
g[r/l{t)]. At distances r <C l{t) from a typical reference 
point inside the cluster the correlation function should 
obey the Porod law: g{^) = 1 — fc^ with a constant k 
of order unity. At /(t) < r < L{t) g(^) ^ ^^"^ (see 
Rcf. Q), a power-law tail with the same exponent as in 
C{r,t^O). 

In order to check these predictions, we solved Eq. (|^) 
numerically on a domain 2048 x 2048 with no-flux bound- 
ary conditions. The accuracy of the numerical scheme 
was monitored by checking the (approximate) conserva- 
tion law (^ which was found to hold with an accuracy 
better than 0.2% for t > 3. 

We used 10 different DLA clusters js^] as the initial 
conditions. These clusters (like the one shown in Fig. 1, 
upper left) had a radius of order 10'^. To prevent frag- 
mentation at an early stage of coarsening, the clusters 
were reinforced by an addition of peripheral sites, similar 
to Ref. 1^ . The average fractal dimension of the initial 
clusters, determined from the averaged pair correlation 
function, was 1.75. 

Introducing the density p(r,t) — (l/2)[u(r,i) -|- 1], we 
identifled the cluster as the locus where p{r,t) > 1/2. 
Typical snapshots of the coarsening process are shown in 
Fig. 1. One can see that larger features of the FC grow 
at the expense of smaller ones. At late times the cluster 
shrinks, in agreement with the prediction of DSI [ p2[ . 
A similar shrinking is evident in the pictures obtained 
in Monte Carlo simulations of area-preserving interface- 
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controlled coarsening [Q , although the authors of Ref . |^] 
did not comment on it. 

To characterize the coarsening process, several quanti- 
ties were sampled and averaged over the 10 initial condi- 
tions: 

1. The cluster area. 

2. The (circularly averaged) correlation function, nor- 
malized at r=0: 

r( i\ < p{r' + Y,t)p{r',t) > 
= <pHr\t)> • 

3. The coarsening length scale l{t), computed from 
the equation C{l,t) = 1/2. l{t) can be interpreted 
as the typical width of the cluster branches. 

4. The cluster perimeter A{t) computed by a standard 
algorithm Q. 






FIG. 1. Evolution of a DLA cluster undergoing an in- 
terface-controlled coarsening in a globally conserved system. 
The upper row corresponds to t = (left) and 12.6 (right), 
the lower row to t = 126.4 (left) and 1856.6 (right). 

The cluster area was found to be constant with an 
accuracy better than 0.5% for t > 10, and better than 
0.15% for t > 100. Hence, area preserving motion by 
curvature, Eq. (||), provides an accurate description to 
this regime. Figure 2 shows that, at late times {t > 100), 
C{r,t) approaches a scaled form. The scaled function 
has a long-range power-law tail with an exponent D — 
2 (the same as in the initial condition), see the inset 
of Fig. 2. Noticeable is the absence of any additional 
dynamic length scales, in a striking contrast to the locally 
conserved fractal coarsening jll]. The dynamics of l(t) 
is shown in Fig. 3. The same figure shows a pure t^^^ 
power-law line (serving as a reference for the expected 
late-time behavior) and a corrected power-law fit l(t) = 
lo + br with a = 0.49, b = 1.2 and = 5.0. 




FIG. 2. Scaling form of the correlation function C(r, t) for 
time moments t = 400.0, 587.0, 1264.8 and 1856.6. The inset 
shows the same data on a log-log plot. The solid line, serv- 
ing as a reference, represents a power-law with an exponent 
D-2 = -0.25. 
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FIG. 3. The coarsening length l(t) versus time (circles). 
The solid line is a corrected power-law fit: l{t) = lo + bt" with 
a — 0.49, lo — 5.0 and b = 1.2. The dotted line represents 
a pure t^^^ power law. The inset shows the time-dependent 
effective dynamic exponent ao{t) versus l/l{t). The solid line 
is a linear fit. 

The inset of Fig. 3 illustrates a different method |Q| of 
determining the dynamic exponent. One defines a (time- 
dependent) effective exponent: ao{t) = d\nl{t)/d\nt. 
Under the normal scaling assumption, one can deter- 
mine the "true" dynamic exponent by plotting ao (t) ver- 
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sus l~^{t) and extrapolating it to i ^ cx), that is to 
l~^{t) — > 0. The numerical value of ao{t) is computed 
from ao{t) = log^g [l{10t)/l{t)]. This procedure yields 
a = 0.50. Therefore, l{t) exhibits normal scaling. 

The same procedures were used for an analysis of the 
dynamic behavior of the cluster perimeter A(i). We 
found that A~^{t) exhibits the same normal scaling: 
A~^(t) ^ i^/^. Irisawa et al. 0| reported an exponent 
0.38 for the inverse perimeter. Their graph shows, how- 
ever, that the effective exponent increases at late times. 
We believe that a careful analysis of their data would also 
lead to an exponent of 1/2. 

Thus, all predictions following from the DSI hypoth- 
esis: the normal scaling of l{t), shrinking of the FC 
and scaling behavior of correlation function (including 
its power-law tail), are confirmed by numerical simu- 
lations. We therefore conclude that globally-conserved 
interface-controlled fractal coarsening exhibits DSI and 
normal scaling. This behavior stands in contrast to 
the breakdown of scale invariance observed in diffusion- 
controlled coarsening of FCs , where the order pa- 
rameter is conserved locally. The mechanism of scal- 
ing violations in locally-conserved systems is not entirely 
clear at present, therefore a comparison between the two 
types of systems can be instructive. We relate the dif- 
ference in scaling behavior to an important (and simple) 
difference in the character of transport. Global trans- 
port, characteristic for interface-controlled systems, is 
uninhibited by Laplacian screening effects typical for lo- 
cally conserved systems. Therefore, large-scale dynam- 
ics is always present in globally-conserved systems, but 
is suppressed in locally-conserved ones. This difference 
is observed already in a very simpler setting of an area- 
preserving coarsening (shrinking) of long slender bars. In 
the locally-conserved case the bar acquires a dumbbell 
shape, while its initial width remains (almost) constant 
and represents a relevant length scale until late times 
Iprf . On the contrary, in the globally conserved case the 
shrinking bar has a finger-like shape, and its dimensions 
are changing on the same time scale [^sf . 
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Sciences and Humanities. 
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